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Abstract 

Important questions concern the existence of excitonic strings in organic 
compounds and their signatures in the photophysics of these systems. A 
model in terms of Hard Core Bosons is proposed to study this problem in 
one dimension. Mainly the cases with two and three particles are studied for 
finite and infinite lattices, where analytical results are accessible. It is shown 
that if bi-excitonic states exist, three-excitonic and even, n-excitonic strings, 
at least in a certain range of parameters, will exist. Moreover, the behaviour 
of the transitions from one exciton to the biexciton is fully clarified. The 
results are in agreement with exact finite cluster diagonalizations of several 
model Hamiltonians. 

PACS Numbers: 71.10.-w, 71.20.Rv, 71.35.-y 
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I. INTRODUCTION 



Excitonic molecules - or biexcitons - are known and well studied since several decades 
in conventional semiconductors. Their formations due to the coupling of a pair of exciton 
and their photophysics properties have been well established both theoretically and exper- 
imentally [I]]. The outstanding facts are related to the giant oscillator strength associated 
with the one exciton-biexciton transition; we may cite, for instance, the giant two photon 
absorption, shown theoretically by Hanamura and observed first on CuCl by Gale and 
Mysyrowicz There, five or six orders of magnitude more than for typical two-photon 
interband transitions are observed. 

In organic compounds, the interest for excitonic bound states - or excitonic complexes 
or last, to use a terminology coming from the Bethe Ansatz technique, excitonic strings - is 
rather new since the main studies started about five or six years ago. They concern three 
types of compounds: some organic Charge Transfer solids (CT), J or H-aggregates (J-ffA) 
and Conjugated Polymers (CP). 

The organic CT solids considered here, consist of planar aromatic Donor (D) and Ac- 
ceptor (A) molecules alternately arranged along a one- dimensional stack. A large distance 
between molecules of the same stack (more than 3 A) gives a rather small hopping integral 
compared to the characteristic Coulomb repulsion. Moreover, a very large distance between 
neighbouring chains ensures the quasi-one dimensionality of these systems. A prototype 
is given by the anthracene-PMDA (pyromellitic acid dianhydride); for this compound, by 
comparison between theoretical calculations on small clusters and experimental results of 
differential transmission spectroscopies, the existence of biexcitons - and even triexcitons 
- has been demonstrated recently 0]. This result stays, for the moment, the only clear 
evidence for excitonic strings in organic compounds. 

J and H-aggregates are also, in many cases, one dimensional stacks of organic molecules. 
But, on the contrary to the previous class of compounds, they have no ionic character 
and the attractive exciton-exciton interaction is smaller. They are typical organic solids in 
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the sense that their low-lying excitations can be understood in terms of Frenkel excitons. 
These compounds are claimed to be the ideal ones to get biexcitons but, without convincing 
demonstration for the moment - possibly be due to a too small binding energy 

Conjugated polymers are very complex compounds. As far as we are concerned with 
low-lying excitations, we may consider them, at first approximation, as an one dimensional 
interacting electron gas, which strongly interacts with the lattice too [0]. Moreover, for a 
more realistic description, in many cases, three dimensional effects and inter- and intra-chain 
disorder should be considered. Because of these difficulties, despite a large amount of work, 
even the nature of the primary excitations remains subject of intense controversies [§. Any- 
way, some Photoinduced Absorption (PA) experiments P,p!0| and Two Photon Absorption 
(TPA) experiments [p!T] , p!2| show some features which are possibly due to biexcitons. 

Only a few theoretical investigations of excitonic strings in organic compounds have been 
done 0-§|,|12|-0. Most of them are based on exact diagonalizations of small clusters using 
different kinds of model Hamiltonians depending on the system under consideration P-|6|,p!3 



From these results and based on physical reasoning, it has been argued that, when biexcitonic 
states exist, the intensity of the optical one exciton-biexciton transition may decrease when 
the binding energy of the biexciton increases; moreover, this intensity is shown to be less 
important than the one of the one exciton-two excitons (two free excitons) transitions and, 
even, is supposed to be independent of the system size for sufficiently large binding energies 
0. Therefore, if the tendencies observed on finite clusters remain the same for the infinite 
system, on the contrary to the situation observed in some inorganic semiconductors , it 
may be difficult to create biexcitons by TPA. A more appropriated experimental technique 
to observe excitonic strings would then be the Differential Transmission Spectroscopy used 
by the way in Throughout this manuscript, we often refer to these numerical works 
for comparisons, especially concerning the intensity of the transition from one exciton to 
biexciton states. 

In this paper, we propose a model for studying excitonic complexes in organic com- 
pounds. Based on physical arguments, we think, this model is relevant for the three classes 



of compounds enumerated above. For organic crystals, as the J and H aggregates, our model 
is equivalent, in the low energy sector, to an effective XXZ spin one-half Heisenberg model 
already used to study this problem p|,p!^-p!6[] . For the CT solids and CP, some extensions 



are necessary because the proper excitons are not of Frenkel type. They are realized by 
introducing effective Hard Core Boson particles with a finite extension on a lattice, which 
mimics the exciton extension. Next, we propose an effective Hamiltonian in terms of these 
particles which captures the essential physical ingredients to describe excitonic strings in 
organic compounds. This model is considerably easier to study than the traditionally em- 
ployed models and therefore some analytical investigations are possible. In this paper, we 
study the proposed model for two and three particles in a finite and infinite lattice. Our 
main interests concern the behaviours of the binding energies of the excitonic strings with 
the parameters of the model and of the one exciton-biexciton transition oscillator strengths; 
are they giant as in inorganic semiconductors or, on the contrary, small as suggested by exact 
diagonalizations on small clusters 0,0? This question is crucial in order to clarify the role 



of hypothetical biexcitons - or even more extended excitonic strings - in the photophysics of 
organic compounds. 

The paper is organized as follow. In section II, we introduce our Hard Core Boson 
Hamiltonian. In section HI, IV and V, we study this Hamiltonian for one, two and three 
particles respectively. 

II. MODEL 

First, for each class of systems mentioned in this study, we emphasize the main charac- 
teristics which are important for our purpose. 

J and H aggregates are typical Organic Crystals. The ground state of these systems is 
the tensorial product of the molecular ground states. The low-lying excitations are Frenkel- 
type excitons where the hole and the electron of a monoexcitation are located on the same 
molecular site PJT^- p!B| , P^ . We can draw a picture of these excitations in the following way: 



4 



if we represent symbolically the ground state of such a quasi-one dimensional system as 
...MMMMMM..., where M stands for a Molecule in its ground state, the Frenkel excitons 
are the Bloch states of local configurations such as ...MMM^ MMM... where the super- 
scripts — and + are for the electron and the hole, respectively; is then a monoexcited 
molecule. 

Organic CT solids are stacks of Donor (D) and Acceptor (A) molecules. We already 
mentioned that they are, with a good approximation, quasi-one dimensional systems with 
a nearest-neighbour hopping term small compared to the Coulombic term. The strong 
coupling regime is then a good starting point to consider such systems. Moreover, if one 
consider only a short range Coulombic term, their ground state would be a Charge Density 
Wave depicted roughly by the unique configuration ...DAD AD A... where D molecules get 
two electrons and A molecules none. The low-lying excitations would then be the excitons 
where an electron is transfered from a D to an A molecule; they would be Bloch states of 
excited configurations such as ...DAD^A~DA... (more details could be found in ref 0]). We 
get a so-called Charge Transfer Exciton (CTE). For a more realistic Coulomb potential, the 
picture proposed above will remain reasonable, at first approximation. 

Conjugated Polymers are quite different from the other systems just considered since, 
obviously, they are not molecular crystals. On the contrary, they are usually described as an 
one dimensional electron gas governed by the Pariser-Parr-Pople (PPP) Hamiltonian which 
includes both long range Coulomb interaction and, semi-classically, electron-phonon interac- 
tion 0. Numerical studies of this Hamiltonian give, for a proper choice of parameters, exci- 
tonic states of small radius, typically a few monomeric units [p!2| , p!7| . On the other hand, the 
ground state may be analyzed in terms of basis sets completely localized on the monomers. 
Then, the ground state of finite clusters is shown to contain mainly intermonomer charge 



fiuctuations of very short range - especially the nearest-neighbours ones [|20[. According to 
these specificities, an effective model has been proposed very recently pT| , p2| . Starting from 
the PPP Hamiltonian, this model reaches a kind of molecular-crystal description for the 
low-lying excitations in conjugated polymers. In this new approach the vacuum (or ground 



state) is still not the simple molecular-type of ground state, proper to the Simpson related 



models p3| , p^ , but a matrix-product like state composed by local configurations extended 
at most, over two monomers. Excitons are then composed mainly by local configurations 
where either the electron and the hole are localized on the same monomer or distant by 
only one monomer; the weight of the configurations where the electron and the hole are 
more separated decreases exponentially in the wave functions. Therefore, the radius of the 
exciton, r, defined as an average distance between the hole and the electron is rather small. 
If we assume r = 1, the following pictorial representation can be done: if we represent the 
ground state as ...MMMMMM..., where M is a Monomer, the excitons will be the Bloch 
states of such configurations {...MMM+ M' MM...) ± {...MMM' M+ MM...). Because of 
the electron-hole symmetry, the charge transfer of the electron on the right side and on the 
left side must be combined either symmetrically or antisymmetrically. In contrast to the 
two previous cases, excitons have here, a structure. 

In conclusion, in these three cases, excitons are Bloch states of very localized excitations. 
The vacuum and the nature of these local excitations vary, of course, following the system 
under consideration: crystal molecular ground state and Frenkel excitons for J and H 
aggregates, CDW and CTE for CT solids and finally, matrix-product state and some kind of 
CTE but including the two possible arrangements of the charges, for CP. However, all these 
very localized excitations can be thought in terms of effective Hard Core Bosons (HCB) 
extended over a few bonds - with an extension of the order of the exciton's radius. We 
adopt such a view in this work, and describe the dynamics of these excitations within the 
same model which we introduce now. 

In this study, we consider only one dimensional lattices with N sites, a site being a 
molecular-site in the case of molecular crystals H] and a monomer-site in the case of CP 



21| , |22[] . As we already mentioned, this is a good approximation for a large class of compounds 
studied here. However, the generalization to higher dimensions should be possible and the 
subject of further studies. 

Excitons are extended over r sites, the radius of the excitons, which varies from zero to 



a few units in the cases considered here. In this study, to mimic the excitons, we introduce 
HCB particles extended over one bond (r = 1). This 'minimal' extension is adopted more 
by a pictorial preoccupation but, also, in accordance with the work developed in ref pl| , p2| . 
Indeed, at the thermodynamic limit where we will derive all the analytical formulae, this 
extension plays obviously no role, since it remains small compared to the system size. More- 
over, for finite size cluster calculations, the qualitative behaviours shown in this paper could 
also be obtained by choosing r = 0, r = 2 or even other 'small' r. The restriction to r = 1 
does not affect the conclusions but, if necessary, more (or less) extended particles could also 
be considered without any difficulties. Working with extended particles, instead of particles 
localized on a site, does not imply technical complications. 

The Coulomb interaction between excitons must be included in order to describe excitonic 
strings . For that purpose, the relative position between the electron and the hole, the two 
elementary constituents of an exciton, becomes important. We introduce then two species 
of HCB extended over one bond: 

• The so-called Right-bosons (R) which mimic nearest neighbour electron-hole pairs with 
the electron on the right (see figure |I|.a). 

• The so-called Left-bosons (L) which also mimic nearest neighbour electron-hole pairs 
but with the electron on the left (see figure ^b). 

The R and L particles obey the following Hamiltonian 

H = LJo EniRiRn + LiL^ ) - J EniRiRn+l + h-C.) - J Y.n{LiL^^^ + h.C.) 

-aEniRiL^ + h.C.) -VEniN^ - N^)iN^+2 - ^^^a) 

RI^^ and L^j^) are the destruction (creation) operators of the R and L bosons at site n, jV^*-^\ 
the operators number of HCB R (L) at site n. ujq is the excitonic energy [ujq > 0). The 
J term is the hopping term of the HCB (J > 0); it is directly related to the exciton band 
width given by 4J. a is an effective term for virtual interactions which couple locally R 
and L bosons (a > 0). The R and L particles are effective representations of more complex 
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particles, invoking a more sophisticated level of description; these "high level" particles 
possess a structure which can be very complicated as it is the case for most of the CP. 
The interaction a reflects the existence of virtual processes exchanging the relative position 
of the electron and the hole of an exciton; this can be due, for instance, to kinetic terms. 
Indeed, by a second order process involving a n.n. hopping term, it is readily possible to 
exchange the relative position of the two building particles of an exciton. The last term of 
(H), y (y > 0), is the interaction between excitons restricted to nearest-neighbours. We 
have to distinguish between two situations: for two excitons of the same species, RR or LL, 
the nearest charges are of opposite sign so that the resulting interaction must be attractive; 
the situation is reversed for two excitons of opposite species, RL or Li?, where the resulting 
interaction is repulsive. 

For structureless excitons, Frenkel or CT-excitons, the definition of two species of degen- 
erate HCB makes no sense. To study strings of such excitons with the general Hamiltonian 
(0), it is sufficient to: (i) identify the R and L particles and (ii) to take the limit a = 0. 
To be more precise, R and L particles may exist in the case of CT solids, but with a large 
energy difference which results, in our model, in a very small a; therefore, it is a good 
approximation to simply neglect this term in this case. 

In the following, we will study the Hamiltonian (|l]) for one, two and three particles, 
which are the cases of experimental interest. 



Excitons appear as Bloch states of symmetric or antisymmetric combinations of R and 
L bosons, corresponding to two different symmetry classes 



where | > is the state without any boson. The separation in energy between this two 
states is given by 2a, the (+) state being lower. 



III. ONE EXCITON STATES 




1 



(2) 
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In the case of J (H) aggregates and CT solids, we identify R and L bosons and get only 
one excitonic state 

|-,A;>=|0> , I +, A; >=| A: >= , ^ V e^'^'^Rj I > (3) 

In all cases, exciton states result in a very intense peak in the linear absorption spectrum 
for k = 0, while the electron-hole continuum is almost not visible. This very unusual feature, 
if we think about conventional semiconductors, is an important characteristic of these one 
dimensional compounds |plip2| , p5| . This is easily understandable since the radius of the 



exciton is very small and the ground state contains mainly short range charge fluctuations 



21,22 



Now, we write down the elementary transition moments due to a localized HCB, rhcp, 
in the case of CP and rhoc, in the case of J (ff) aggregates and CT solids (oc is for Organic 
Crystal); these quantities will serve in the following to express the transition moments 
between one exciton and two exciton states 

= ^ < I ef(i?t + Lt ) I > 



--< I efRi I > 



where e is the charge of the electron and r the position operator. 

The intensity of the ground state-one exciton transition at A; = is proportional to 
I fhoc P or I fhcp P, depending on the system under consideration [^. In the case 

of CP, only the excitonic state | +,0 > is observable in linear absorption. However, the 

state I — , > could become important, as it is the case in electroabsorption experiments for 

instance. 

IV. TWO EXCITON STATES 

The wave functions with two excitons are written as 

I ^2 >= J2 '^cna2ini,n2) \ ni,n2 > (5) 

ni<n2—l 
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where ai and cf2 are R or L bosons, | ni,n2 > the ket relative to the situation where 
the particle cxi is in ni and (T2 in n2 with ni < 77.2 — 1, ^Po-io-al^i' ^2) is the corresponding 
amplitude. 

As usual, we treat separately the movement of the center of mass and of the relative 
position of the two particles. Moreover, we have to specify the nature, R or L, of the 
particles located on rii and n2, it is then convenient to use matrix formalism. We write 



^2>= E 



JQ{ni+n2) \ 



\n I ni,n2 > 



(6) 



n>l,ni 



with n = n2 — ni (n > 1) and 



(7) 



where the small left (right) indices refer to the nature of the left (right) particle, R or L 
bosons. 

We apply Periodic Boundary Conditions (PBC) to the wave functions which give the 
following constraints 



Q = ^ , j = '^,-;N {j an integer) 
[v]^n+N = l^]n where n = n2 - Ui 



(8) 



Moreover, we work with Hard Core Bosons extended over one link so that we must 



preserve 



]o = Ml = [0] 



(9) 



where [0] is the zeroth four component vector. 

With the wave function (H), the eigenvalue equation is written 



E[ip]n = -2JcosQ[/]{[v?]„-i + M„+i} - am[ip]n - <2V"[A]Mr 



(10) 



10 



where [/] is the 4 by 4 identity matrix, and 



^0 10 1^ 



10 10 
10 1 
10 10 



[A] 



V 



-1 



(11) 



With the equation (pnD, we have reached an impurity hke problem where free particles, 
belonging to four different continuum, interact with impurities pinned at n = 2. First, by 



applying a local unitary transformation, we solve (10) away from the impurities 



mn = [u]H 



(po{n) 



( ^ 

V2 
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V 2 



iPrl{n) 

ipir{n) 



\ 



We can now rewrite equation (^) in terms of we find 



(12) 



E\ 



-2JcosQ[/]{[(^]„_i + [ip]n+i} - 2a[a2Mn - Sn,2V{[aiMn + [(TsMn} (13) 



where 



1 

1 



(14) 



Without the impurity like states at tt, = 2, there is no mixing between lpq, ip^, ip+ and 
V9_ components; with the scattering of particles at n = 2, only the and components 
are mixed. 
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From equations (^) and (|T2D, it is clear that only the if+{n) and the '^-{n) components 
of [0\n are relevant to study the optical responses of two exciton states. Since we are mainly 
interested in spectroscopic properties of biexcitons, we consider only these two components 
in the following. The (po{n) and ^Qin) components have no physical meaning for structureless 
excitons, but could be important for CP; their characteristics could be calculated, following 
exactly the same way as the one presented here. 

We write, following H. Bethe p6| . 



ip+{n) = A+e'^" + S+e"'^" and ip^{n) = A„e^" + B.e"^" (15) 

which give the following expressions for the energy where J{Q) = 2 J cos Q and where we 
choose 2u!q as reference 

E = -2a~ J{Q) {e" + e"^) and E = 2a - J{Q) (e'' + e"") (16) 



y4_|_, A_, and B_, fi and rj are the constants to be determined. Note that, if in (p!5D 
/i and r] are purely imaginary numbers, we get free excitons. On the contrary, if they are 
purely real numbers, we get bound-states, the biexcitons. 

The ansatz (^) gives obviously the exact solution of (|I^) without impurities (V^ = 0); 



with impurities, ( |l5l) must fulfill the following constraints in order to remain solution of (|13| 



J{Q)V^{1)-Vv42)=0 ^^^^ 

Moreover, the PBC imply 

A+ = e'^^B+ and A. = e""^ (18) 

In the following, we give the solutions of the above mentioned problem for different 
choices of parameters. 
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A. a = 0, = 



In this trivial case, the two excitons are free; we get two degenerated continuum for ip+{n) 
and V5-(n) within the energy range [—2J{Q),2J{Q)]. It is easy to show that a very intense 
linear absorption takes place for Q = and /i = ?7 = only; then the transition moment of 
the one exciton - two free exciton transition at the edge of the continuum behaves as 

I M^i exciton—2 free excitons \ ~ \ ^^cp/oc \ (-^*^) 

plus corrections in For a 7^ and/or V 7^ 0, this strong absorption remains but slightly 
shifted toward the high energies. 

This transition is similar to the ground state - one exciton transition; indeed, the excita- 
tion energies and the behaviour of the transition moments are the same in these two cases. 
This characteristic was already used in [^,^ to explain the lack of bleaching signal in PA 
experiments for Charge Transfer and Frenkel exciton systems. 



B. Q = 0, y / 

This case concerns CT solids and J or H aggregates for instance. Without the a term, our 
model is equivalent to a standard interacting Frenkel exciton Hamiltonian already studied 
numerically and analytically |T^-|TB[. Our results, obtained with a different technique, 
are in accordance with the conclusions reach by these various studies. 

Since y 7^ 0, bound states - or biexcitons - may exist; in this subsection we consider the 
biexcitonic states only and then and rj of equations ([T5|) , as real quantities. 

When a = 0, equations (|l^) imply that fi = rj. From equations together with the 



PBC ([ISD , we get at the thermodynamic limit, the following set of equations 

(20) 

A2 _ = 

where A = 

13 



Additionally to the continuum, if A < 1, i.e V > J{Q) = 2 J cos Q, we get with the 
second equation of (|20| ) and equations (p!6|), two bound states with energies given by the 
following simple expressions 

By using a two-particle Green function approach |T^, these equations have already been 



obtained for Q = in ||T5[ and for every Q in [|1J,|1^ together with the biexciton condition 
V > J{Q). 

The first equation gives a bound state below the continuum of two free exciton states, 
this is the biexciton of interest for us; the second equation gives another bound state on top 
of the continuum. In the following, we consider only the first state and calculate its wave 
function and the transition moment associated with the one exciton-biexciton transition. 

First, we write A± = e'^/^ and S± = e^'^. With the PBC (|TB|) and equations (^, we 
get the expression for 5 

-5 = N\n\ (22) 
Then, with equations (|T5|) and (pOD, we get 



=7\Acosh(-lnA[n - — ]) (23) 



valid in the interval l<n<A— 1. A/" is the normalization constant given by 



(24) 



The first equation of (pOD shows that the wave function of the lowest bound state is a 
symmetric mixture of the tp^ and components resulting in configurations with two like 
particles RR or LL. On the contrary, the highest bound state contains only configurations 
with two unlike particles RL or LR. We recall that for this case R and L particles are 
identified so that the distinction between the case with RR or LL bosons and RL or LR 
bosons is then a way to treat the case of attractive and repulsive interaction between bosons 
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at the same time. In case of attractive interaction (two bosons of same species), we get 
a bound state below the two exciton continuum; the situation is reversed in the case of 
repulsive interaction (two bosons of opposite species), where a bound state appears above 
the two exciton continuum 0. 

For large systems, the transition moment of the one exciton-biexciton transition is non- 
zero only for Q = and k = {Q, the center of mass of the biexciton, k, the exciton 
momentum); its expression is then given with the help of the elementary transition moment 
(i) by (more precisely, the square of the transition moment is shown here) 

I M.C r= 2 I rUocY^Mn) \'= 2 | ^ " A^/^"^))' (25) 

Going to the thermodynamic limit, we found the simple result 

\Moc\^=A\moc\^^^ (26) 



which is exactly the result found in [p^pi)| by using a two particle Green function approach. 

We recall that the bound states (biexcitons) exist only if < A < 1 and that the 
biexciton binding energy increases when A approaches 0. With this formula, the behaviour 
of the absorption from the one exciton to biexciton state appears very clearly. First les 
us consider its behavior when A approaches or 1, at the thermodynamic limit (equation 
(p6|)): for A ^ 0, the transition moment reaches saturation to 4 | rrToc P which is easily 



understandable since, in this case, the two excitons are strongly bounded together in two 
n.n. sites; for A ^ 1, the transition moment diverges following the asymptotic behaviour in 
4 I riioc P -jzx- Next, we can make some conclusions by studying the equation (pSj). 

• At large and fixed A^, with decreasing A (which means, with increasing the biexciton 
binding energy), the intensity of the transition decreases. 

• At fixed A and sufficiently large A^, Moc is almost independent of the system size. On 
the contrary, the intensity of the transition between one exciton and two free excitons 
increases linearly with the system size (cf. eq. |I^). However, as expected, when A 
approaches 1, | Moc P recovers progressively a linear dependence in A^. 
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Both of these behaviours depict inverse mechanisms than the ones theoretically predicted 



for the ground state-one exciton transition ||21| , p2| , where the intensity of the excitonic peak 
is proportional to the size of the system and increases with the binding energy. 

These behaviours are illustrated in figure |^ where the equation p5| ) is evaluated for 
= 100. There, we see that the intensity of the one exciton-biexciton transition decreases 
dramatically when A decreases (which means, when the binding energy increases) reaching 
saturation for small A. Then, this intensity is independent of which gives in PA and 
TPA spectroscopy an intensity proportional to instead of A^ for the transition to the 
continuum of two free excitons; this last behaviour can be seen by comparison with equation 
(p6|) evaluated at the large A^ limit (dashed line). Consequently, at the thermodynamic limit, 
biexcitons may be observable by spectroscopy experiments, but with the need of very clean 
compounds, otherwise we may expect at first analyse, that the intense transition to two free 
excitons would not permit to detect possible biexcitons. Within our model, no two-photon 
giant resonance can be expected for J, H aggregates and organic Charge Transfer Solids. 
The very same conclusions were already pointed out in ref from numerical calculations 
and in [11^,0 from analytical studies. 



C. a / 0, y / 

This case concerns CP. As in the previous subsection, we look for bound states only so 
that fi and rj are again assumed to be real quantities. The question we want to address more 
specifically here is how the a term will affect the binding energy, the wave function and the 
transition moment related to the biexcitons. 

For a ^ 0, from equations (|T^ and ([T^), we get at the thermodynamic limit the following 
set of equations 

B_ = Ae-'^e^^'S, 

(27) 

A^e** -6-/^ = 

From equations (|T6D and the second equation of (^7]), we get 
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_ ..... 1 



/5 + J/32-(l-A2)(l--) (28) 



1 - AM V ' A2 



where A = and /3 = This expression together with equations (|T6|), gives the energy 
of the lowest biexciton. 

The binding energy of the biexciton, is then given by the following equation 

E,(a) = J(Q)[(e-^ + e'^)-2] (29) 

The study of this expression shows that the biexciton binding energy decreases when a 
increases (see figure ^). This behaviour is not surprising since the a term mixed together R 
and L particles; by doing so, configurations with two unlike particles, RL or Li?, appear in 
the wave function which increase its energy and then reduce the binding energy with respect 
to the case for a = 0. The critical value Oc for which the binding energy becomes zero is 
given by the solution of the equation + = 2, 

The critical value, Oc, increases when A decreases, to diverge for A = 0: for a very large 
binding energy, a has no effect on the bound states anymore. Another interesting quantity 
along this critical line, is the critical value Vc 

- = = ' (31) 

From this equation, we can see without any surprise, than - the critical value above 
which one gets biexcitonic states - continuously increases when [3 increases. The attraction 
between excitons must be stronger when the mixing between R and L bosons becomes more 
important. 

We analyze now the wave function and the resulting transition moment associated with 
the one exciton-biexciton transition at Q = 0. We proceed as in the previous case; we write 
Aj^ = a+e~~, = a+e~ and v4_ = a_e~~, = a_e~ and get with equations 



and (13) 
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V^.i., = AT, (^cosHf,[n - ^]) + ^_iY.,,.,^i,+inX)(N-3) + 2 In A] - ^])) (32) 

where again 1 < n < N—2. Afa is the normahzation constant which could be easily calculated 
if necessary. The first term of (|3^) comes from the ipj^ component, the second term from the 

component. For a = 0, fi = — ln(A) and the previous result is recovered. For a 0, the 
weight of the component decreases smoothly in the wave function. Consequently, the a 
priori unfavourable configurations RL or LR appear in the biexciton wave function, on the 
contrary to the previous case without the a interaction. 

The transition moment of the one exciton-biexciton transition at Q = can be easily 
calculated with the expression of the wave function (^2l), but contains no new information 
compared to what was already pointed out for a = 0. The results obtained in the previous 
subsection remain valid in this case. The only difference comes from the fact that the 
component is non active in linear absorption. Consequently, for the same binding energy, 
if we assume | m^p | = | ?7\,c |, the biexcitonic peak would be a bit more intense for J (H) 
aggregates or CT solids than for CP. For instance, for a = 0, we get | M^p p= ^^""^ ; for 
a 7^ 0, the Lp+ becomes more important than the ip- component in the wave function: the 
intensity of the one exciton-biexciton transition then increases. 

Now, let us take reasonable parameters for Conjugated Polymers. Extracting from ex- 
perimental data (we take roughly the binding energy of the IB^ exciton state observed in 
Polydiacetylene compounds), we get 4 J = 0.5eV By calculating the Coulomb interac- 
tion between two HCB using the Mataga potential (for instance) and for a reasonable 
lattice constant (a = l.sA), we get V = 0.5eV. The a term may be evaluated from the 
energy difference between the IB^ exciton and the nA~ "ionic" exciton. These states are 
represented here, in an effective way, by | +,0 > and | — ,0 > respectively (see eq. 0). This 
quantity could be estimated, for instance, from exact resolutions of the PPP Hamiltonian 
for small clusters or, in a better way, from electroabsorption experiments. There, the state 
I +,0 > shows a red-shift with no change of the spectral line and a quadratic dependence 
on the applied field. This quadratic Stark effect of the exciton can be reasonably well re- 
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produced with the state | — ,0 > only, neglecting, in particular, the continuum effects p5 . 
Within this assumption, it is then very easy to estimate a. We take here a reasonable value, 
a = 0.2eV. With this choice of parameters, we get a binding energy of 0.12eV for the 



biexcitons which is in accordance with experimental results PJTl[]. By the way, it is easy 
within our model to get such a value for the binding energy for other values of a or /and J 
by changing V, which is more flexible. 

In figure ^, a calculation for a finite cluster of ten sites is shown. The parameters are 
A = 3.33 and (3 = 0.2 which give a typical spectrum. BE is for the peak associated with the 
BiExciton and, FE, for the peak associated with the more intense two Free Exciton state 
(even if such concepts are not well defined for finite clusters) . The BE peak is a bit smaller 



than the FE one, which is in accordance with the results of ref. |jT^ from exact calculations 
using the Extended Hubbard Hamiltonian, a short version of the PPP Hamiltonian. At 
the thermodynamics limit, only the FE peaks will survive, slightly shifted toward the low 
energies. 

To conclude with the two particle case, we may summarize our main results. For a = 0, 
biexcitons exist if < A < 1; the transition moment of the one exciton-biexciton transition 
shows inverse behaviour compared to the ground-state-one exciton transition: it decreases 
when the biexciton binding energy increases and, for sufficiently large A^, is independent of 
the system size. For a ^ 0, the conclusions aforementioned remain valid at the condition that 
a < ac for J and V fixed. Above this critical value, biexcitons fall in the continuum. Finite 
size cluster calculations give qualitatively similar results as other calculations performed 
with more complicated Hamiltonians. 

V. THREE (AND MORE) EXCITON STATES 

For more than two particles, the question about the integrability of the Hamiltonian (|1]) 
becomes important. To answer this question, it is sufficient to analyze carefully the case 
with only two particles P7| , p5| . For that purpose, let us first rewrite the Hamiltonian (|I|) in 
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a more convenient way by introducing two new operators 

(33) 



Then, the Hamiltonian becomes diagonal in a 

H = {uj^~a)Y.nBlB^ +{uj^ + a)Y.nAiA^ - JE„(St5„+i + /i.e.) ^^^^ 

To illustrate this new picture, note that the states with one exciton (Q) are the Bloch 
functions of one B or one A particle. 

Let us now consider the case with two like particles {BB or AA) located in ni and n2 
(ni < 722 — 1). Because the last term of ([3^ ) exchanges the color [B or ^4) of the two 
particles, the corresponding amplitude in the wave function is written as 

^(m, ns) = C(5, ki] B, A;2)e*^i"i+*^2n2 + c(B, k^, B, kAe'^^''^+'^^''^ + 

(35) 

C{A, k[; A, ]^^)eik[m+ik'^n2 ^ ^^ J^^yik'^n^+ik[n2 

the constants C being complex numbers determined in the previous section by solving the 
Schrodinger equation. 

The system is invariant by translation, hence the total momentum is conserved: ki + k2 = 
k'i + k'2 = Q. We then write ki = Q — q, k2 = Q + q and k[ = Q — k, k'2 = Q + k. 

Now, let us do a scattering 'experiment'. Starting, for instance, with two B particles, 
after the scattering processes we get two A particles. The collision being elastic, the energy 



(p!6|) must be conserved 

-2a- J{Q) (e« + e"^) = 2a - J{Q) {e'' + e''^) (36) 
Depending on the a value, we can make some conclusions. 



If a ^ 0, q k (see eq. ^ for instance). Hence, there is a change in the momenta 
during the scattering processes. In other words, diffractive processes are involved 
during the scattering which is sufficient to conclude that the model defined by (|I]) (or 



4])) is non-integrable 
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li a = 0, q = k. There is no diffractive processes anymore; then, as we look for the 
lowest bound-states, we will see that the model (1) may be expressed, in that case, as 



an effective spin-1/2 XXZ Hamiltonian which is well known to be integrable p7|-|30 
Rigorously, to conclude about the integrability of this model in the general case, the 
study of the two particles S matrix has to be done We will not attempt such 

study here. 

In the following two subsections, we will study both cases, starting with the simplest one 
where a = 0. 

A. a = 

For a = 0, by interacting with each other, particles do not exchange their color; the 
interaction is reduced to an exchange of momenta alone. In such case, it is better to consider 
the Hamiltonian (|1]) expressed in terms of R and L particles. 

Since R and L particles are HCB, an n-electron configuration is subdivided into n sepa- 
rated regions. Each particle stays in its own area of the chain; within optical terminology, 
only refiections are allowed. The Hamiltonian is studied independently for each configura- 
tion of color which corresponds to a configuration of pair-interactions, attractive or repulsive, 
depending on the nature of the neighbouring particles. In the general case, the problem is 
very complex. However, the interesting case for us is given by the configurations where all 
the particles have the same color. In this case, the model can be written in the form of a 
XXZ ferromagnetic Hamiltonian for spin-1/2 which has been solved exactly by using the 



Bethe-Ansatz p|J2|-pO| 



First, our goal is to show the equivalence of our model (|1|) in the low energy sector and 
the XXZ model. For simplicity, we consider here particles without extension (r = 0). Of 
course, it is possible to work with our defined R and L particles, but the comparison with the 
XXZ model would then claimed some additional notations unnecessary for our purpose. By 
the way, in the second subsection of this section, we will show more formally the equivalence 
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between the two models working with the extended particles. 

Second, since we consider the simplest case with n-particles like - for instance, the states 
with n R particles - the Hamiltonian of the system may be expressed in terms of spin-one- 
half operators: S"^ = +| being for a site n without any particle and = — |, for a site 
occupied by a particle (an exciton). The model becomes then equivalent to the spin-1/2 
XX Z model with an 'attractive' anisotropic term P| 

H = 2ujo^ - 2J^ -(S'+5'^+i + - ¥^^{3^ - - -) (37) 

n n n 

The Ground state of ( P7| ) is the state with only up-spins, \GS >= \... |-[-|-[-| ... >. 5"- creates 
a down-spin (a R particle or exciton) localized at site n, S:^\GS >= \... TTTinTTT ••• >• We 



introduce here, as usual, the parameter of anisotropy A = jj p8|-p0| 



This model has been subject of intense studies since several decades, and many results 
are known about it |]28|-pO|]. Concerning the bound states, or more precisely, the so-called 
"string-states", their properties depend on the anisotropic parameter A. Within the "string 
hypothesis", the following results were found [p^ -|3T||. 



If A > 1, the bound states of n-spins, here n-excitons, exist without any restriction 
on n [p8H 3l|| . The energy of a n-excitonic string is then given by the following formula 



first appeared in 



=nu;o-nV-J sinh 7 '"'^^^"^ (38) 

smh nj 

with cosh 7 = A and Q being the wave number of the n-string center of mass. For 
n = 2, the energy ( |2TD is of course recovered. 

If < A < 1, the n-excitonic strings may exist but with a strong restriction on n 

{n - 1)0 < TT (39) 



with cosG = A, < < |. Moreover, these excitonic strings fall in the continuum 
spectrum. 
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We have seen in the previous section, that biexcitons exist if A > 1. In this case, from 
the above conclusions about the XXZ-model, we can say that n-excitonic strings will also 
exist and that their binding energies will increase with n (|38D; this result is in agreement 
with numerical studies on small clusters 0]. However, as we have shown in detail for the 
biexcitonic case, it is not easy to observe such states experimentally, since the particles are 
tightly bounded inside the complex, resulting in a lack of oscillator strength. For < A < 1, 
excitonic strings may still exist but with an energy falling in the continuum. In this limit, 
the excitons are less bounded and, therefore, it could be interesting to study the signatures 
of these states under the influence of an electric field. We leave this question to further 
work. 



B. a / 

For a 7^ 0, the interactions between particles involve both exchange of momenta and 
exchange of color. By a simple argument, showing the appearance of diffractive terms in 



the two particles scattering case, the non-integrability of the model p^ ) has been shown. 
Therefore, to study n-excitonic strings with n > 2 one needs to do some approximations. 
In the following, we propose a trial wave function for the three particle case build from the 
Bethe Ansatz solution of the XXZ-Heisenberg model within the string hypothesis. 

Let us consider the expression ( P^ of our model expressed in terms of A and B particles. 
For the three particle case, eight configurations have to be distinguished depending on the 
color of the particles and on their relative positions; we denote these configurations as criO"2a"3, 
where (Tj is for the color [A or B). These configurations are separated into two disconnected 
channels: on one hand, we have AAA, BBA, BAB and ABB; on the other hand, BBB, 
AAB, ABA and BAA. We consider only the latter channel which is the lowest in energy. 

The wave function for the triexciton is expressed as a linear combination 

I i/jt >= ao I ^0 > +ai I i'l > +0-0 \ '^c > +ar \ A > (40) 
where x = 0,r,c,l stands for the BBB, AAB, ABA and BAA states respectively (r, c, / for 

23 



right, center and left refer to the position of the B particle in the three particles complex). 
ttx are the coefficients to determine. Every component | > is assumed to have the Bethe 
ansatz structure 



*x >= J2 i'x{ni,n2,n3) \ ni,n2,n3 >^ 

l<ni<n2<n3<N—l 



(41) 



where, more precisely, the particles are again separated by one bond, at least, and with 



(42) 



The kf are the quasi-momenta and the 6'^- the phase terms pq^O|. The summation runs 
over the 3! permutations of the indices. We apply periodic boundary conditions to the wave 
functions, '>pxij^2, n^jUi + N) = -02.(^1, n2, ns), and we get the well known equations ^6 



ik: 



(43) 



Since we are looking for excitonic complexes, we are interested in complex solutions of 
these equations. Following Ovchinnikov we assume that the only relevant phase terms 



are 6*23, then 



Nkf = 2^Af + 



Nkl = 2ti\% - ^fg 



(44) 



where are integers between and — 1. It follows 



Nlrakf = Im^fg 

Nlmkl = Im(^f3 - Of^) => 



Nlmk^ 



2Im6'f2 = A^(ImA;f — ImA;f ' 



2Im^f3 = iV(Im/cf - Imk^) 



(45) 



Next, we assume 6^2 = Im^f2 < ^23 ~ Ii^^23 < without loss of generality; then, 
at the thermodynamic limit, among the summation (^21) only one term remains with an 
exponential accuracy in 



V'x(rii,n2,n3) = e*'=?"i+*'=2'^2+ifcfn3+e-2+ef3 



(46) 
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We may notice, that the used assumption about the phase terms 6'f-, is rigorously correct in 



the case without a 31 



The wave function ( |40| ) with the simphfication (^Uf), must satisfy the following constraints 
due to the V term, in order to be an eigenfunction of (0): 

• for 712 = ^1 + 2, (Vna) 

aQj{^Q{nx + + 2,713) + ^0(^1, ''ii + 1,^3)) - aiV'^i{nx,nx + 2,713) 
aiJ{^i{nx + + 2,713) + ^z(rii,7ii + 1,713)) - ao^*o(?^i, + 2,713) 
acJ(\l'c(ni + 1, Til + 2, 713) + \l'c(ni, rii + 1, 713)) - ar'V'^rijix.nx + 2, 713) 
J(\E'r-(m + 1, rii + 2, 713) + \l/r(?ii, iii + 1, 713)) - a^V'^cirix, rii + 2, 713) 



: 

= 

= 



(47) 



for ?i3 = 7i2 + 2, (Vrii) 

ao J(^o("'l, ""-2 + 1, ?l2 + 2) + V&ol'^il, ""-2, n2 + 1)) - ar^^r(?T-l, ""-2, 7l2 + 2) = 

arJ{^r{nx,n2 + l,7i2 + 2) + '^r{nx,n2,n2 + 1)) - ao^^ol^T-i, "-2, ''T-2 + 2) = 

acJ(^c(m, ?^2 + 1, 712 + 2) + ^c(ni, 712, ?l2 + 1)) " 0/^ (ill , 7l2 , 7l2 + 2) = 
aiJ(^^/ (771, 772 + 1, 772 + 2) + */(ril, r72, 712 + 1)) " ac'^^^c(ril, r72, 712 + 2) = 



(48) 



Moreover, by symmetry, we have k\ = k\, 6'[^j+i = 9\^^i,\fi and 
previous set of equations together with the expression (^) gives 



ai. Then, the 





e''°2 






+1 




ik^ 
e 3 






+1 


- Vf^^'^Xj 


■ik^ 

e 2 






+1 




e 3 






+1 



(49) 



to which we add the momentum conservation law 



(50) 



Q being the momentum associated with the motion of the center of mass of the excitonic 
complex. This system of equations is obviously non-soluble in the general case. Hence, we 
add an additional simplification to the trial wave function. 
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If a = 0, /cf = and therefore = |,Vx; it is easy to verify that the system (^91) 
becomes equivalent, in that case, to the corresponding equations for the XXZ Heisenberg 
model |31|| as we already pointed out in the previous subsection. Then the system (^Q]) 
leads to the energy (pH]), for A > 1. Within the string hypothesis, and k^ are complex 
conjugate, k2, a purely real quantity. 



Enlightened by the exact results for a = we assume for a 7^ 

k{ = k^^ + e 

kl = kl (51) 

U — UO _ f: 

where e, a real number, is given by the following equation obtained from the energy conser- 
vation law 

- 3« - ^(e'^° + e-""') =a- f](e''=' + e-""'') (52) 

i=l i=l 

The left- hand- side is the energy of the BBB component, the right- hand- side, the energy of 
the three other components. 

With this assumption on the momenta fcf, the solutions stay very close to the Bethe 
Ansatz form; the binding between the two bordered particles of the complex (with momenta 
ki and /cf), ensures the cohesion of the three particle bound states. The set of equations 



composed by the two first equations of ( ^9[ ) and equations (|50D, ( [5T| ) and ( |52D are then solved 
numerically. 

It is important to realize that our approximation is variational. This important statement 
could be checked in the following way. One starts with the trial wave function given by (^Op 
where the components | '^x > are assumed to be of the form defined by the equations 
(1^) and (^) together with the restrictions (|50|) and (|5^) . The coefficients of (j^), and 
their complex conjugates, are determined variationally by minimizing the functional 

J^({a^, <}) =< I I > < I > (53) 

where —E is the Lagrange parameter for the normalization constraint of the wave function. 
We get 
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, , , , , , 

7— = ^ ao < ^0 hJ'o > +cir <'^o \ H \ -^r> +ai <'^o \ H \ '^i>= aoE < > 
da*Q 

(54) 

and similar equations for all the other minimizations. At the end of the day, we obtain 
nothing else than the usual secular equations and one way to solve them, following A. Bethe 
2^, is to impose the constraints given by the sets of equations (^T]) and (^S]). Our procedure 



is then variational and the results found are an upper bound for the problem. 

Typical results are shown in figure (^). The binding energy of the biexciton got in the 
previous section and the one of the triexciton, obtained with our variational calculation, are 
represented for Q = and A = 3.33 as a function of /5 = j. The triexciton binding energy 
is defined with respect to the continuum of one biexciton plus one free exciton which is the 
lowest in energy. The following remarks can be made, (i) The triexciton binding energy for 
a = is the exact one. (ii) The binding energy of the triexciton is larger than the one of 
the biexciton. (iii) A critical value of a exists also for the triexciton; it is larger than ac ( |30D 
for the biexciton. It follows that triexcitons may exist without biexcitons (in the sense that 
Ef, > 0). We may suspect the same behaviour for the binding energy of n-string for any n: 
the binding energy and the critical a may be larger for a n-string than for a (n-l)-string in 
accordance with the results for the XXZ model, (iv) Since our result is variational, we get 
an upper bound for the triexciton binding energy so that our conclusions are qualitatively 
correct. 

To illustrate the behaviour of the binding energy of a n-string at strong a, we study 
this limit perturbatively. For a >> V, the Hilbert space of ( |5^ is naturally separated in 
two subspaces: the first, lower in energy, with configurations without A particles and, the 
second, higher in energy, with configurations with at least one A particle. The first one is 
associated with the projector Pq, the second one with the projector (Fq + -Pr^ = !)• The 
Hamiltonian may be written as 

H = PoHPo + P^HP^ + PoHP^ + P^HPo (55) 
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By a Schrieffer-Wolf type of canonical transformation p2| , |33| , we want to derive an 
effective Hamiltonian H = PqHPq + Pr^HPr, + W, with P^HPr, = ^,W being the effective 
term to be determined. The study of the n-excitonic strings will then be performed by 
diagonalizing the projection of this model into the space without A particle, Heff = PqHPq. 
Introducing the unitary transformation U = (with S*^ = — S"), and following the method 



described in detail in ||33|, we found the simple effective Hamiltonian valid in the strong 
coupling limit where a >> V, at second order in ^ 

H^ff = {uo - «) E BiB^ - 2JE i(^i^n+i + h.c.) -^Y. BIB^bU.B,,^, (56) 

n n ^ n 

We have recovered, once again, the spin-1/2 XXZ Hamiltonian (|37D. However, the 
anisotropic parameter is now given by A = instead of and the interpretation of 
the spin operators is different: 5"^ = +| is still for a site without any particle but, = — |, 
is now for a site occupied by a S particle instead of a i? particle. 

In the effective model, the attractive interaction is reduced by the a term and, for 
the strong limit (a >> V) where this derivation is valid, one reaches a XX-Heisenberg 



Hamiltonian which is well known to be equivalent to a spinless free fermion system 



Hence, in this limit, the excitonic strings do not exist which is consistent with the results of 
our variational calculation for a > ac- The states with free B particles are then preferred. 



VI. CONCLUSION 

In this paper, we propose a simplified model to describe excitonic strings in quasi-one di- 
mensional organic compounds including organic Charge Transfer solids, J- and if-aggregates 
and Conjugated Polymers. For all these compounds, the excitonic states are characterized 
by small radius. Hence, our model starts with the definition of Hard Core Boson particles 
extended over one bond which describe, in an effective way, the excitons in one dimensional 
organic compounds. 

There are two possibilities in one dimension to place the electron and the hole, the two 
elementary constituents of an exciton. Either the electron is on the right or on the left of the 
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hole; we introduce two kinds of HCB which illustrate these two situations: the Right-hosons 
and the Left-bosons. 

The proposed model in terms of these effective particles contains four terms: the excitonic 
energy (uq) , which separate in energy the space with n particles from the space with n — 1 
particles; the kinetic energy of the HCB (J); the local interaction between the two species 
of bosons (a), which illustrates the fact that there exist effective interactions which can 
exchange the relative position of the hole and the electron of an exciton; and last, the 
interaction between two excitons which are attractive between particles of the same species 
and repulsive between particles of opposite species (interaction of intensity V). The kinetic 
energy and the a-interaction combined with the repulsive two body interaction act against 
the attractive two body interaction - they tend to delocalize the excitons on the contrary 
to the two particle attraction which tends to create bound states. The cases with two and 
three particles are mainly studied in this work. We summarize briefly our results considering 
the case where Q, the momentum of the center of mass of the excitonic complexes, is zero. 

The case with two particles is solved exactly. For a = 0, biexcitonic states exist if ^ > 1. 
Then, the intensity of the transition from the one exciton state decreases dramatically with 
the increase of the binding energy, to saturate and become independent of the system size. 
This behaviour makes the biexcitonic state difficult to observe at the thermodynamic limit. 
For a 7^ 0, a critical value etc is found, above which biexcitonic states do not exist. Otherwise, 
the behaviour of the oscillator strength remains qualitatively unchanged. 

For a = 0, the n particle case can be solved in the low energy sector. Then, the model 
is equivalent to the spin one half XXZ Heisenberg model which is exactly solvable in one 
dimension. For ^ > 1, it is well known that n-excitonic strings exist with a binding energy 



which increases with n [p8|-pl|]. For a ^ 0, the model is non-integrable. The three particle 
case is then studied using a variational ansatz, which may be extended to a more general 
case with n particles. Again, a 3-string is found with a binding energy larger than the one for 
the biexciton. We also found a critical value for a, larger than the ac valid for the 2-strings. 
From numerical studies on small clusters 0], it has been argued that n-strings may exist 
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for any value of n in organic compounds. Our work gives some confidence to this statement. 
Indeed, with our HCB Hamiltonian, if biexcitons exist, triexcitons will exist with a larger 
binding energy; moreover, if a = 0, n-strings will also exist with binding energies which 
increase with n and, we believe, from our present results, that it will also be the case for 
a ^ 0. 

The next important question concerning the n-string peaks in PA and TPA experiments, 
is fully clarified for the biexciton (the most important case) and confirms some conclusions 
of previous works [^,|6i|T3|-[l6[| . The intensity of the transition from the one exciton to the 



biexciton decreases with the binding energy and, for sufficient binding energy, becomes 
independent of the system size. Hence, the PA intensity to the biexciton is proportional to 
^/N instead of for the two free exciton case. This property comes from the fact that the 
two excitons are tightly bounded within the complex, rendering the 2-string states difficult 
to observe for infinite systems (clean samples are needed) and more accessible experimentally 
for small oligomers P,p^|. 



In the case of Conjugated Polymers, we think some additional studies are needed going in 
two different directions. First, in the context of PA experiments, excitations extended over 
two molecules have been invoked P,p5|; extension of our work to two dimensions or, at least. 



for two coupled polymers would be suitable. Second, Conjugated Polymers are disordered 
systems with several possible sources of disorder which result, in practice, in the vague 
definition of the so-called conjugation length [^. As we have already mentioned, since the 
transitions to biexcitons are much less intense than the ones involving two free excitons, we 
may expect at first analyse - and it was implicitly our point of view all along this manuscript - 
that a sufficient disorder will render the biexcitons not observable. However, some relatively 
recent results about the two interacting particle problem in a disordered medium show 
very unexpected behaviours which may have dramatic consequences for our problem, and 
possibly in a way reverse to our intuition. In the light of the intriguing results of ||3^, the 
effects of disorder on transition moments should be clarified for a correct interpretation of 
PA and TPA results. We leave these considerations to further works. 
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FIG. 1. Pictures for a) Right-bosons and b) Left-bosons. The arrows show the position of the 
electron, the hole being on the other side 
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FIG. 2. Intensity of the one exciton-biexciton transition for N=100 (full line) and at the ther- 
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modynamic limit (dashed line) as a function of A = ^ {a = 0, Q = 0) 
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FIG. 3. Normalized transition moment for the one exciton - 2 exciton transitions for a finite 
cluster (N=10), A = 3.33 and /3 = 0.2. BE (Bound Excitons) is the biexcitonic peak; FE (Free 
Excitons) is the more intense peak due to two free excitons. At the thermodynamic limit, only FE 
will survive, slightly shifted toward the low energies. 




FIG. 4. Negative of the binding energies for the bi- and tri-excitons for A = 3.33 in function 
of/3=f (Q = 0) 
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